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ABSTRACT 

The use of Pontryagin's Maximum Principle for the large-angle slewing of large 
flexible structures usually results in the so-called two-point boundary-value problem 
(TPBVP), in which many requirements (c.g., minimum time, small flexible amplitude, and 
limited control powers, etc.) must be satisFicd simultaneously. The successful solution of 
this problem depends largely on the use of an efficient numerical computational algorithm. 
There are many candidate algorithms available for this problem (c.g., quasilinearization, 
gradient, and shooting, etc.) In this paper, a proposed algorithm, which combines the 
quasilinearization method with a lime shortening technique and a shooting method, is 
applied to the minimum-time, three-dimensional, and large-angle maneuver of flexible 
spacecraft, particularly the orbiting Spacecraft Control Laboratory Experiment (SCOLE) 
configuration. 

Theoretically, the nonlinear TPBVP can be solved only through the shooting method 
to find the "exact" switching times for the bang-bang controls. However, computationally, 
a suitable guess for the missing initial costates is crucial because the convergence range 
of the unknown initial costates is usually narrow, especially for systems with high 
dimensions and when a multi-bang-bang control strategy is needed. On the other hand, 
the problems of near minimum time attitude maneuver of general rigid spacecraft and fast 
slewing of flexible spacecraft have been examined by the authors through a numerical 
approach based on the quasilinearization algorithm with a time shortening technique. 
Computational results have demonstrated its broad convergence range and insensitivity 
to initial costate choices. 

Consequently, a combined approach is naturally suggested here to solve the minimum 
time slewing problem. That is, in the computational process, the quasilincarization method 
is used first to obtain a near minimum time solution. Then, the acquired converged initial 
costates from the quasilincarization approach arc transformed (tailored) to and used as 
the initial costate guess for starting the shooting method. Finally, the shooting method 
takes over the remaining calculations until the minimum-time solution converges. The 
nonlinear equations of motion of the SCOLE arc formulated by using Lagrange 7 s 
equations, with the mast modeled as a continuous beam subject to three-dimensional 
deformations. The numerical results will be presented and some related computational 
issues will also be discussed. 

’ Research partially supported by NASA Grant NSC.-I4I4 and supplements. 
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INTRODUCTION 


* Future space missions (target acquisition, target tracking, 
and surveying multiple targets, etc.) require. 

• large-angle rotational (attitude) maneuver (slew) ; 

• 3 -dimensional (3-D, 3-axis) maneuver; 

• large flexible spacecraft maneuver; 

• minimum time maneuver . 

* Application of Pontryagin’s Maximum Principle to the nonlinear 
slewing problem: 

I. Non-Minimum-Time Slews: 

1. 3-D Rigid Spacecraft . 

Junkins, Turner, Vadali, Wie, Bainum and Li, etc. 

2 . 2— D (Single— Axis Rotation) Flexible Spacecraft 
Turner, Junkins, Vadali, Chun, Thompson, Bainum and Li, 

etc. 

3. 3-D Flexible Spacecraft (SCOLE) 

Bainum, Li and Tan. 


IX. Minimum-Time (Near-Minimum-Time) Slews: 

1. 3-D Rigid Spacecraft 

Bainum and Li, Vadali, Wie, etc. 

2. 2-D Flexible Spacecraft 

Singh, Junkins, Vadali, Byers, Bainum and Li. 

3. 3-D Minimum-Time Flexible Spacecraft; Using 
Quasilinearization Method and Shooting Method : 
present paper. 
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OUTLINE 


1. 3-D Dynamics of Flexible Spacecraft 

• State Equations 

2 . Time Optimal Control Problem Formulation 

• Two-Point Boundary-Value Problem (TPBVP) 

3. Quasilinearization Method for Near Minimum Time Slew 

4. Shooting Method 

5. Initial Costate Transformation 

• Scale Factors 

• Combined Algorithm 

6. Numerical Examples 

7. Conclusions 
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3-D DYNAMICS OF FLEXIBLE SPACECRAFT 


I. 3-D Deformations of the SCOLE Mast (modal superposition) : 

U='L^ i (z) a i (t) , v=Ei] i (z) a i (t) , <J> = EC i (z) a i ( t) (l) 

i i i 

where 

U , V — bending in x and y directions; 

0 - torsion in z direction; 

£ i , and Ci - modal shape function vector components; 

- a scaled modal amplitude associated with the ith mode; 
z - coordinate. 

II. State Equations : 


Q = 
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where 



(A +B a ) 0 ) + (Cp)to + Da + (E +F a ) u 


(4) 


q is the 4 x l quaternion vector. 


B a = [B ltt 

! ■®2 a i • • • ! ] , 


C R = [C X B 

\ C 2 & | C 3 6], F a = 

! ■^ 2 ^ !•••! *^*9®] t 

A - , , C L , 

D, E, F i - constant 

matrices; 

u = t f lx 

■^ly ^ lz i ^ 2x f 2y i 

■^3x ^3y i *^4x 
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TIME OPTIMAL CONTROL PROBLEM FORMULATION 
Two-Point Boundary -Value Problem (TPBVP) 

Initial States and Final Requir ed States: 

g(0) , a (0) , a) (0) , P (0) ( 5 ) 

g( t f ) , a( t £ ) , w(t f ), p(t f ) ( 6 ) 

Cost Function : 

t f = /"'(licit (V 

Jo 

Saturation-Bounded Controls : 

| u i | ^ u ib , i = 1, 2, 9. ( 8 ) 


Hamiltonian : 

H - 1 + Y t P + ^ T [ U +B a ) O + Cp&) +£>« + {E +F a ) u] (9) 

p t y f 1=[1 1 X 2 ] T “ costate vectors associated with q, a, <o, B. 
costate Equations (by Pontryagin's Maximum Principle): 

(10) 


v . -M . - P-X - - (f. t D U (11) 

da 

= --§£ = -4 idp - +s a )]w - (c p ) t i ( 12 > 

da> 2 

i 2 = -Jjj[ = - Y - (CpA) o) < 13 > 


Constraint Condition ( 
Optimal control : 


a terminal condition to determine 
H = 0, o <. t Z t f 


sign[ (E+F a ) T A] it i = 1, 9 . 


t f ) : 


(14) 


(15) 
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QUASILINEARIZATION AND TIME SHORTENING METHOD (QTS) 

The near minimum time slewing problem can be solved by using 
an iteration approach based on the quasilinearization method. 
Different Cost Function (Quadratic Cost Functional) : 

J = [ tf (u T Ru)dt (16) 

2 Jo 

where R are weighting matrices, t f is the slewing time. 

Same Costate Equations Eqs. (10)-(13). 

Different Control Expression : 

Unconstrained Optimal Control : 


§H 

du 


0 , 


u = -R' 1 (E +FJ T X 


(17) 


Constrained Optimal Control : 



U iC 

if 

\ u ic\ < u ib 


u i = 1 

u ib sgn(u ic ), 

if 

l U icl * U ib 

(18) 

= ~[R 

'He +F a ) T A] i , 

i = 

1, 2, 9. 

(19) 


t f can be obtained by sequentially shortening the slewing time. 
Motivation : 

Is this bang-bang control the same as that obtained by using 
the shooting method? (Do these controls have the same time 
histories?) If the answer is yes, the results from the QTS 
approach may be used as the starting solution for the shooting 
method. (Here, we use the numerical results to prove the 
equivalence. ) 


299 



SHOOTING METHOD 


Formulation of the TPBVP : 

X(t) = F[X(t) , u( t) ] , 0 <; t s t f (20) 

X(0) = K[D] < 21) 

L [X( t f ) , D, t f ] = 0 (22) 

u i = -sign {g i [X(t)] }, i = 1, m- ( 23 ) 


D - the n x 1 unknown initial costate vector; 

L[X(t f ) , D , t f ] - (n+1) x 1 terminal constraint vector; 
g L (i=i , m) - the switching functions. 


Tnitial Boundary Conditions Cor rection Process: 

To satisfy: L[X(t f ) , D , t f ] = 0, D and t f need to be corrected at 
each iteration: 


D ' 

tl 
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dL 
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dD 
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L[X ik) ( t f ) , D (k) , tj* 1 ] 

D ( k ) and - the values of D and t f at iteration k; 

scalar a k (0 < < 1) is chosen as: 

pi [D (k) , tj W ]| 


a* = min 


1, 


|| [ 6D (Jc) , 6tj k) ]« 


>, 0<P<1 


(24) 
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(26) 


(27) 
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A general algorithm suggested by Lastman has been used in our 
calculations and can be recast into the following block diagram: 



Summary : 

I. Difficulty in initialization for the present nonlinear, multi- 
input system control problem. Improper choice of D and t f will 
result in singular correction matrix, and program diverges. 


II. Advantages and disadvantages of the two methods: 

• Quasilinearization method and time shortening technique 
has good convergence properties; 

• Shooting method generates more accurate final results, 
but is sensitive to starting solution; 

• A combined technique is needed. 
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INITIAL COSTATE TRANSFORMATION 


Different Costate Solutions : 

Although the costate equations are the same, the initial 
costates from these two different formulations of the same problem 
are different: 

1. Initial costates, p( 0) , y( 0), and 1(0) from the QTS 

method are large; 

2. p(o) , y ( 0) , and 1(0) from the shooting method are finite. 
Assumed Relation Between Them : 

p( 0) * iqp(0) , 1,(0) * W0) , 


Y (0) * k 2 y(0) , 1 2 (0) = ic 2 l 2 (0) . 


(28) 


where k 1 < 1, and k 2 < 1 are scale factors to be determined. 


Scale Factors 

Assumed Eigen-Axis Rotation for Rigid Spacecraft : 

to = e 0 , o) = e 0 (29) 

e = [ej^ e 2 e 3 ] T ~ a unit vector representing the eigen-axis, 

0 - the rotation angle about this axis. 

Resulting Four Related Equations (from rigid dynamic equations) : 

e 0 = f 0 2 + I _1 Bu ( 3 °) 

/ - 3 x l constant vector; and 

0 = e T f 0 2 + e T I~ 1 Bu (31) 

• Let "p" = the "principal" axis among the axes 1, 2, and 3, 
about which the rotation requires the largest t f ; 

• Let "4" = Eq. (31). 
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Further Simplification of Equations : 


6 = b i v 1 , b i > 0 , | v i | s 1 , i = 1 , 2 , 3 , 4 . 

"Average” Values of the Initial Costates : 

p a = 2/ (b a Q f ) V*. \ a = l/b a , t fa = 2 ( e f /b a ) V* 


b a = k p b p + k t b it k p + k t = l 

0 f - the required rotation angle about the eigen-axis. 


/c 1 and k 2 : 

*i = P«/|P p (0) |/ k 2 = \ a /\\ p (0) | 

The Initial Costates for Starting the Shooting Method : 


p(0) 

= k x p ( 0 ) , 

Xi(0) = k 1 X 1 (0) , 

Y(0) 

= * 2 y(0) , 

X 2 (0) = k 2 k 2 ( 0) . 
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NUMERICAL EXAMPLES 


Given Slewing Conditions (for all cases considered here) : 

1. Rest-to-rest slews, i.e., 

0 ) (0) =0, to ( t f ) = 0, 

a (0) =0 , p(0)=0, a(t f ) =0, p<t f ) =0. 

2. Three (3) control variables are used. 


Example 1 (a scaled rigid spacecraft, A = : 

g(t f )=[ .877582561, .434965534, .142572492, . 142572492 ] T , 

g(0) = [ 1 . 0 , 0.0, 0.0, 0.0] T , 0 f =1 rad, I = Diag(1.0, 0.9, 0.6), 
R = Diag (1.0, 0.7, 0.4), p = 1, k x = k p = k A = 0.5. 

OTS Method Results : 

t f = 1.8 sec. By transformation, the initial costates for 
starting the shooting method are obtained: 



-1.67968 ' 


-.705333 

p(0> - 

-.248420 

, A (0) = 

-.0955727 , 


.415782 


.0481507 


t f =1. 71209 
(s) 


Shooting Method Results : 

p = o.l, solutions are obtained in 6 iterations (to 5 digits) : 



-1.74008 ' 


-.770403 

* 

o 

II 

-.267243 

, A* (0) = 

-.115614 


.462349 


. 0606796 


t;=l. 76403 
(s) 


The converged values of the switching times are: 


u i 

u 3 

u 2 

U 1 

u 3 

u 2 

t i(s) 

0.314356 

0.701830 

0.874531 

1.18114 

1.53158 
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Example 2 ; 


Maneuver of the rigidized SCOLE model. The scaled inertial 
matrix (set J 33 = 1.0): 


I = 


.16902 -.001061798 .01619427 

-.001061798 .9948471 -.007354633 

.01619427 -.007354633 1.0 


• Simultaneous 75°, 30°, 45° slew about roll, pitch, yaw axes; 

• u = three torquers on the Shuttle, u ib = 10,000 ft-lb; 

• R = Diag ( 1 , 1, 1) ; 

• k p =k 3 =0.75, k 4 = 0.25. 


Average Values : 


p a =1.01466, A a =.461915 
Initial Costates for Starting Shooting Method : 



.00238330 ‘ 


-.00709927 

p( 0) = 

1.24510 

, A(0) = 

.456237 


-1.01466 J 


-.461519 


The t f = 1.6407 sec from the QTS method is used as t£ 0) in the 
shooting method, p = 0.0035. The final converged initial costates: 



.00237568 ' 


' -.00773848 ' 

P*( 0) = 

1 . 32922 

, r<o) = 

. 496600 


-1 . 08389 


-.502920 


=1.64066 

(s) 
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Time (sec) 



Time (sec) 



Time (sec) 



Time (sec) 


Fig. 2. Rigid SCOLE Maneuver (Scaled Model) 
75-30-45 (deg) (Shooting Method) 


306 







Example 3 : 


The maneuver of the both rigidized and flexible orbiting SCOLE 
model is considered (using the original SCOLE challenge 
parameters) . 

• 90° slew about roll axis. Although the expected motion is 
"single— axis" rotation, the minimum— time dynamic maneuver 
process is not necessarily a single-axis rotation because of 
the offset inertia distribution of the SCOLE model. (The three 
axes for the three control torquers are not principal axes) . 
Therefore, the present slew is a 3-D slew. 

• u = three torquers on the Shuttle, u ib = 10,000 ft-lb; 

• R = Diag ( 1 . E-4 , l.E-4, l.E-4); 

• k p =k l =1 * 

• Two flexible modes (the first and the second) are included. 
Average Values of Initial Costates : 

p a =24.7475, A a =120.242 
Initial Costates for Starting Shooting Method : 

7 ELC0J ylOj 1 1 I0J i 2 101_ 

(Rigid) (Flexible) (Rigid) (Flexible) 

•00000E0 .00000E0 -.90512E-3 .12024E3 .12024E3 .10717E-1 

. 24747E2 .24747E2 .46390E-2 .67830E1 .51236E1 -.35993E-1 

-.19219E0 - . 11478E0 .11671E2 .10969E2 

. 25248E1 .23758F1 

The t f = 27.3992 seconds from the QTS method is used as t^ 0) in the 

shooting method. The final converged initial costates: 

7 eioj um A^IO) i 2 X0j_ 

(Rigid) (Flexible) (Rigid) (Flexible) 

. 00000E0 . 00000E0 -.87221E-3 .11894E3 .11587E3 .10331E-1 

. 24422E2 .23767E2 .44563E-2 .67130E1 .49375E1 -.34685E-1 

- . 18880E0 - . 11001E0 .11544E2 .10570E2 

.24825E1 . 22818E1 


307 



The Hamiltonian, H, is observed as a constant during each iteration 
and is iteratively reduced to the final value: 

H= 1 . 2000E-9 


The time histories of the slews are plotted in Fig. 3 to 
Fig. 6. In these figures, the results for the attitude angles, the 
mode amplitudes, the control torques, and the switching functions 
are presented, whenever applicable. 

Fig. 3 shows the rigid SCOLE maneuver by using the QTS method 
and Fig. 4 shows the rigid SCOLE maneuver by using the shooting 
method. The results show that the solutions by using the two 
methods are very close. It is also noted that, during the slew, the 
yawing control, u 3 , switches twice consecutively before other 
controls (rolling control u 1 or pitching control u 2 ) switch. 

Figs. 5 and 6 show the flexible SCOLE maneuver by using 
the QTS and the shooting methods, respectively. Again, the results 
from both methods are close. Due to the inclusion of the flexible 
modes, the switching number for every control is tripled or even 
more (23 for u 3 ) compared with the results for the rigid SCOLE 
maneuver. The modal amplitudes are very small and the associated 
vibration of the reflector of the SCOLE and the "Line of Sight" are 
also very small. 
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(b) Controls (1E4 ft-lb) 



Fig. 3 Rigid SCOLE Maneuver, 

Roll Angle ■ 90° (QTS Method) 
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(a) Attitude Angles (deg) 



6 10 16 20 26 30 

Time (sec) 


(c2) Pitch Controls (1E4 ft-lb) Switching Func. 



"I 1 I I I L 


10 16 20 25 

Time (sec) 


(c3) Yaw Controls (1E4 ft-lb ) S. Func. (IE-4) 

1 
0 
-1 

o 5 10 15 20 25 

Time (sec) 



Fig. 4 Rigid SCOLE Maneuver, 

Roll Angle - 90° (Shooting Method) 
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(a) Attitude Angles (deg) 





Fig. 5 Flexible SCOLE Maneuver, 

Roll Angle ■ 90° (QTS Method) 





(a) Attitude Angles (deg) 



(b) Flexible Modes 



Fig. 6 Flexible SCOLE Maneuver, 

Roll Angle ■ 90° (Shooting Method) 
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(cl) Roll Control (1E4 ft-lb) 


Switching Func. 



Time (sec) 


(c2) Pitch Control (1E4 ft-lb) Switching Func. 



O 5 10 15 20 25 30 

Time (sec) 

(c3) Yaw Control (1E4 ft-lb) Switching Func. 



0 5 10 15 20 25 30 


Time (sec) 

Fig. 6 (Cont’d) Flexible SCOLE Maneuver, 
Roll Angle s 90° (Shooting Method) 
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CONCLUSIONS 


1. The QTS method is stable for relatively coarse choices of the 
unknown initial costates, and the shooting method is not. 

2. The QTS method usually results in very large values of the 
costates which may lead to the numerical overflow in the 
calculation process, if the "exact" (numerically) switching 
times are to be found, while the shooting method does not have 
this problem and exact switching times can be obtained 
iteratively. 

3. A technique to combine these two methods is proposed. 

4. The estimated initial values of the costates, p( 0) , y(°)/ and 
A(0), based on the solution from the QTS method, are very 
close to the converged values of these parameters in the 
shooting method and hence the convergence of the shooting 
method has been improved. 

5. The costates from both methods are proportional. 

6. The control histories from both methods are the same and may 
imply the uniqueness of the control for the slewing problem. 

7. The application of this method to the minimum time maneuver of 
other flexible spacecraft is suggested. 
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